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Quantum interference between energetically close states is theoretically investigated, with the 
state structure being observed via laser spectroscopy. In this work, we focus on hyperfine states of 
selected hydrogenic muonic isotopes, and on how quantum interference affects the measured Lamb 
shift. The process of photon excitation and subsequent photon decay is implemented within the 
framework of nonrelativistic second-order perturbation theory. Due to its experimental interest, 
calculations are performed for muonic hydrogen, deuterium, and helium-3. We restrict our analysis 
to the case of photon scattering by incident linear polarized photons and the polarization of the scat¬ 
tered photons not being observed. We conclude that while quantum interference effects can be safely 
neglected in muonic hydrogen and helium-3, in the case of muonic deuterium there are resonances 
with close proximity, where quantum interference effects can induce shifts up to a few percent of 
the linewidth, assuming a pointlike detector. However, by taking into account the geometry of the 
setup used by the CREMA collaboration, this effect is reduced to less than 0.2% of the linewidth 
in all possible cases, which makes it irrelevant at the present level of accuracy. 

PACS numbers: 32.70.Jz, 36.10.Ee, 32.10.Fn, 32.80.Wr 


I. INTRODUCTION 

Quantum interference (QI) corrections were intro¬ 
duced in the seminal work of Low [T] , where the quantum 
electrodynamics (QED) theory of the natural line pro¬ 
file in atomic physics was formulated. These corrections 
go beyond the resonant approximation and set a limit 
for which a standard Lorentzian line profile can be used 
to describe a resonance. Often referred to as nonreso¬ 
nant corrections US!, they contain the full quantum in¬ 
terference between the main resonant channel and other 
non-resonant channels, which leads to an asymmetry of 
the line profile. Therefore, the fitting of spectroscopy 
data with Lorentzian profiles becomes ambiguous, since 
it leads to energy shifts that depend on the measure¬ 
ment process itself US]. A careful analysis of the lim¬ 
its of the resonance approximation is thus mandatory for 
high-precision optical and microwave spectroscopy exper¬ 
iments. 

The first calculation of QI was made for the Lamb shift 
in hydrogen-like (H-like) ions and for the photon scatter- 
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ing case [T]. It was found to be relatively small, of the 
order of Jqi/F « a{aZ)‘^, compared to the linewidth F. 
Here Jqi is the line shift due to QI, a is the fine structure 
constant and Z is the atomic number. Thus, for some 
time, little interest has been addressed to these correc¬ 
tions in optical measurements of H-like ions. Since the 
late 1990s, high-precision measurements of the Is — 2s 
transition frequency in hydrogen renewed the interest in 
these corrections |S!8|. Numerous theoretical calcula¬ 
tions of QI were then made for H-like ions [mini HO] 
with astrophysical interest m and application to laser- 
dressed atoms [12]. QI has also been studied in other 
atomic systems and processes during the last decades, 
mainly because near and crossed resonances of hyperfine 
states can enhance the QI effects HMH- QI effects dSl- 
IT9| have been shown to be responsible for discrepant mea¬ 
surements of the helium fine structure [SQlEI] and the 
lithium charge radii determined by the isotope shift |S]. 
The lithium experiment |5] gives a beautiful experimen¬ 
tal demonstration of the geometry dependence of the QI 
effect. Very recently, QI effects in two-photon frequency- 
comb spectroscopy have been investigated, too [22] ■ 

In this work, we calculate the QI shifts for 2s^2p tran¬ 
sitions in H-like muonic atoms with hyperfine structure. 
The physical process considered here is the photon scat¬ 
tering of initial 2s states to final Is states, with an in- 
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FIG. 1. (Color online) 2s and 2p level structure (not in scale) for (a) and /i^H (b). Allowed electric dipole 

transitions are also shown. / and F are the nuclear spin and total angular momentum, respectively. 


TABLE I. Calculated energy differences of and 

[231124] states relative to the lowest n = 2 hyperfine 
state, which are 2sf/^°, and 2sf/^^, respectively. r 2 p 

is the respective approximate linewidth. 





GHz 

meV 


2si/2{F = 

1) 

5520 

22.83 


2pi/2{F = 

0) 

51631 

213.53 


2p\/2{F = 

1) 

53521 

221.35 


2p3/2{F = 

1) 

54616 

225.88 


2P3/2(F = 

2) 

55402 

229.12 


Tap [25] 


18.5 

0.0765 


2si/2(F = 

3/2) 

1485 

6.143 


2 pi/2(F = 

1/2) 

49680 

205.46 


2 pi/2(F = 

3/2) 

50182 

207.54 


2P3/2{F = 

3/2) 

52016 

215.12 


2P3/2(F = 

1/2) 

52104 

215.49 


2p3/2{F = 

5/2) 

52286 

216.24 


Tap [25] 


19.5 

0.0806 

p^Re 

2si/2(F = 
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171.40 
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311456 

1288.08 


2 pi/2(F = 

0) 

325633 

1346.71 


2p3/2{F = 

2) 

347860 

1438.63 


2P3/2 {F = 

1) 

353731 

1462.92 


Faf] 


318.7 

0.1318 


Calculated in the present work. 


cident photon energy that is resonant with intermediate 
2p states. Figure recalls the 2s — 2p level structures 
that are considered in this work. The level structures and 
linewidths |23l |25] of the 2p states are given in Table |I] 
The theoretical formalism used here can be traced back 
to recent works EIIMIIIZI- 

Hyperfine states in are separated by several hun¬ 
dred GHz [2^ [28l - [30| . and have linewidths of a few tens of 
GHz [5^. Thus, it is expected that QI plays a small role 
and cannot be responsible for the so-called “proton ra¬ 


dius puzzle”, where a discrepancy of four linewidths was 
observed in the experiments of the Charge Radius Ex¬ 
periment with Muonic Atoms (CREMA) collaboration 
[2I1I22]- Nevertheless, these systematics need to be care¬ 
fully evaluated and quantified, since they have contribu¬ 
tions similar to small QED corrections (e.g. sixth-order 
contributions m) and thus may impact precise deter¬ 
minations of the proton charge radius. Since resonances 
in /x^He+ [331 [M] have been measured recently, values of 
the corresponding QI contributions are presented here, 
too. In the case of /r^H, there is a close energy proximity 
between the states 2p^J'^^'^ and 2p^J^^'^ of « 87 GHz 
and hence, it is expected that QI effects could be much 
higher. 


II. THEORY 


Photon scattering is a two-step process consisting of 
photon excitation with subsequent photon decay, which 
is formally equivalent to Raman anti-Stokes scattering. 
It is described by second-order theories (e.g. Kramers- 
Heisenberg formula [85] . or S-matrix [36]), which overall 
converge to the following scattering amplitude (velocity 
gauge and atomic units) from initial to final states |261 
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where |/), \v), and |/) represent the initial, intermediate 
and final hyperfine states of the muonic atom or ion. 

= E,j — Ei is the transition frequency between \ v) and 
I*). The dipole approximation, ap-e^ (7 = 1 , 2 ) is used, 
where, p is the linear momentum operator and is the 
incident (scattered) photon polarization. The summation 
over the intermediate states \v) runs over all solutions 
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FIG. 2. (Color online) Adopted geometry for photon scat¬ 
tering of incident linear polarized photons with momentum 
fei and polarization ei and scattered photon momentum k 2 , 
uniquely defined by 0. The scattered photon’s polarization is 
not observed in our measurements and is thus not illustrated. 


of the Dirac spectrum of the muonic atom (ion) with 
hyperfine structure. All states are considered with a well- 
defined total atomic angular momentum F, projection 
along the quantified axis m, and total orbiting particle 
angular momentum J [32]; thus the contribution of off- 
diagonal terms (mixing between the and 

states) [35] is considered null. is the full width at half 
maximum for an isolated Lorentzian line of the excited 
state, where we assume T 2 pP =T 2 p. 

The incident photon energies studied in this work com¬ 
prise the near-resonant region of the 2s —>■ 2p transitions. 
This includes only the resonances illustrated in Fig. [l] 
Hence, we restrict the summation v to the 2p states of 
the terms uj 2 s 2 p— 0 Ji [first part of the right side of Eq. ( [l])]. 

Energy conservation leads to Ei — Ef = a ;2 — wi ]26| 
between the initial (Ei) and final {Ef) energy states and 


the energy of the incident (wi) and scattered ( 102 ) pho¬ 
tons; thus only one of the photon energies is independent. 
Using this relation, it is convenient to introduce the en¬ 
ergy sharing parameter u = uii/ujr defined by the fraction 
of the incident photon energy relative to the lowest res¬ 
onant energy uJr of a given muonic atom with initial F) 
(see Fig. [^. 

Motived by the experimental configuration, we con¬ 
sider incident photons having linear polarization and 
non-observation of scattered photon’s polarization, as il¬ 
lustrated in Fig. [^[33]. If we define the scattering plane, 
containing both photon momenta (fci and k 2 ), then a sin¬ 
gle polar angle 6 is sufficient for describing the angular 
distribution of fc 2 - 

The corresponding differential cross section of the am¬ 
plitude in Eq. 0 for all the mentioned approximations 
is given by [33] 


0 y) = _^_ 
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where ui^i = Wi^ijuJr, F^ = T^/uir and Ufi = u + ujfi. 
In Eq. 0 , it is assumed that the initial state of the 
atom is unpolarized and that the level and magnetic sub- 
levels of the final state, as well as the scattered photon’s 
polarization (£ 2 ) remains unobserved in the scattering 
process. are the dipole matrix elements (length 

gauge). Equation ([^ can be further rearranged as a sum 
of Lorentzian components Aj’J''(0,x), and cross terms 

{d,x), similar to Ref. |S]. For our particular ge¬ 
ometry and atomic system, the result is given by 


— (11 fl vi = 

dU ^ (2E, + 1) (a;,, - u)^ + \f,/2)^ 


E 






{^ui ^ 'F y j ^F){fij i/l i U ZT^yl2) 


(3) 


where the second summation over E' and J' runs for 
non-repeated values of Fy and Jy of the first summation. 


The quantities defined by 


mi,Ff ,mf ,Jf ,£2 

and 


2Re 


^i,pf 
"/• U'*2 
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FIG. 3. (Color online) Scattering differential cross-section ([^ at 6^ = 90°and x = 0° for (a) and /r^H (b) versus photon 
frequency. A zoom plot of the two close-lying resonances for is given in the right panel (c). Solid blue lines correspond to 
the full evaluation of Eq. while the dashed red curves account for only the first term in the right side of Eq. § (sum of 
Lorentzian profiles). 


contain all the polarization and geometrical dependen¬ 
cies. {0, X, £ 2 ) and are given in the 

Appendix. 

The differential cross section of Eqs. (§ and (§ con¬ 
tains a coherent summation over resonant excitation 
channels; thus it takes into account channel-interference 
between neighboring resonances. However, as can be ob¬ 
served in Eq. (§ , if cross terms were removed, it re¬ 
duces to an incoherent sum of independent Lorentzian 
profiles. The QI effects are thus included in those cross 
terms. 


III. RESULTS AND DISCUSSION 

In this section, we present and discuss results for the 
QI contribution in several muonic atoms taking into ac¬ 
count Eqs. @ and ffl, first assuming a pointlike detector 
and later the CREMA geometry. The influence of the ge¬ 
ometric and polarization conditions on the QI is well de¬ 
scribed in Ref. [5] and is reproduced in the present work. 
We thus restrict our geometrical settings to the perpen¬ 
dicular observation [9 = 90°) of scattered photons and 
to the case of horizontally and vertically polarized pho¬ 
tons with respect to the scattering plane (x = 0° and 
X = 90°). Figure displays the scattering cross section 
for the 2 s —)■ 2 p —>■ Is processes in and on 

one hand, having the full coherent summation of Eq. ([^ 
(i.e. with QI), on the other hand, having the summa¬ 
tion restricted to only the Lorentzian terms [neglecting 
cross-terms in Eq. ([^]. The peaks correspond to the re¬ 
spective transitions shown in Fig. As it is observed, 
the influence of QI is more noticeable in regions between 
resonances, where no dominant excitation channels exist. 
Close to resonances, the influence of QI is approximately 
equivalent to shifting the peak position, as shown in the 
zoom plot of Fig. 

We determine this shift in each resonance by gen¬ 
erating a pseudo spectrum that follows the theoretical 


profile of Eq. (§ and fitting it with an incoherently 
sum of Lorentzians (as performed in the data analysis 
of the CREMA experiments). Fits are done using the 
ROOT/MINUIT package [40]. All fit parameters (posi¬ 
tion, amplitude, and linewidth) are free fit parameters 
for each transition. The fit range is chosen sufficiently 
large, such that the fit results do not depend on it. The 
shifts of the fitted resonance position, Jqj, normalized to 
r 2 p, for each resonance and muonic atoms are given in 
Table [^ Overall, with the exception of some resonances 
in /r^H, QI produces relative shifts i5qi less than 3% of 
r 2 p. For /r^H and /i^He’*', the shifts are of the order of 
~ 0.2%-1.7% (~36 MHz-310 MHz) and ~ 0.2%-3.2% 
(~ 0.6 GHz—10 GHz) of their linewidths, respectively. 

The observed discrepancy (proton radius puzzle) of 
~ 75 GHz (0.31 meV) jSH |32] at the resonance 

corresponds to 4 linewidths. Hence, it is much larger 
than any possible QI contribution. Moreover, this res¬ 
onance is ~7 times more intense than the closest reso¬ 
nance (see Fig. |^, which minimizes the QI shift 

in this resonance. The values presented in Table [^ for 
^^H and /r^He+ have the same order as the respective 
ones given by the rule-of-thumb for distant resonances 
(^QI °° r|p/4A with A being the energy difference be¬ 
tween two resonances) |16j . Apart from this, relatively 
low intensity resonances, like in ^^He"*", can have 

higher QI contributions due to a high intensity resonance 
nearby. 

F —\ /2 

On the other hand, the resonances ‘^P ^/2 ^ 
^ 7 * 3 / 2 ^^^ in /r^H are more sensitive to QI effects not only 
due to their close proximity (87 GHz), but also due to 
the intensities being comparable within a factor of ^0.7. 
In this case, the QI shifts can be up to 12% of r 2 p 
(- 2 GHz). 

Applying the previously calculated cross sections to 
the geometry of the CREMA setup leads to considerable 
cancellations of the quantum interference effect. 

Figure [^sketches the experimental geometry. Muonic 
atoms (ions) are formed in an elongated gas volume 
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TABLE II. Shift of the line center due to QI contribution 
for /r^He”'" and /r^H. |i) and \v) stand for the initial 

and resonant atomic state, respectively. Values of <5 divided 
by V 2 p are given for a pointlike detector with vertical (dqi, 
X ~ 90°, 6 — 90°) and horizontal (5qi, x = 0°, 0 = 90°) 
polarization cases, as well as for the CREMA setup geometry 



N) 

\y) 

(%) 

(%) 

(%) 


2P3/2 

- 0.8 

1.6 

-0.02 



0.4 

-1.7 

0.01 


^^112 

2.vlT^ 

‘2P3/2 

- 0.2 

0.5 

-0.01 



- 0.6 

1.2 

-0.01 



- 0.5 

0.7 

-0.01 




0.3 

-1.2 

0.01 

p’^He+ 

^^1/2 


- 0.4 

0.7 

0.00 



2P3/2 

0.2 

-0.6 

0.05 


^^112 

‘2ptr2 

- 0.3 

0.5 

-0.01 



- 0.6 

1.4 

-0.01 



‘2pIJ? 

2P3/2 

- 0.1 

0.4 

0.00 



2.2 

-3.2 

-0.02 


^*1/2 


- 0.3 

0.7 

-0.01 




- 0.5 

1.0 

-0.01 




- 3.2 

6.7. 

-0.08 




4.4 

-8.0 

0.11 


^F=ZI2 

^*1/2 


- 0.3 

0.7 

-0.01 



0 F=3/2 

2Pl/2 

- 0.4 

0.9 

-0.01 




- 1.2 

1.9 

-0.03 




- 4.9 

12.3 

-0.13 




0.9 

-2.7 

0.03 


of 5x12x190 mm^ that is illuminated from the side 
using a pulsed laser [HI [42] and a multipass cavity 
[33] ■ The 2p —^ Is photons emitted after laser-induced 
2 s —?► 2p transitions are detected by two x-ray detec¬ 
tors (14x150 mm^ active area each) [33ti3H] placed 8 mm 
above and below the muon beam axis. 

The simplest situation is that the 2s —>■ 2p excitation 
happens in the center of the target (red circle in Fig. |^, 
and the resulting k -2 photon is detected in the center of 
the top detector (black circle). This corresponds to the 
pointlike detector case oi 9 = 90° and x = 90° (i5qj in 
Table |n|. 

However, the 2s —>■ 2p excitation takes place anywhere 
in the muon stop volume, and the photons of the 2p —>■ Is 
decay are detected anywhere on the detectors surfaces. 
We consider here the laser’s propagation (fci) and po¬ 
larization (si) directions being along ei and 63 , respec¬ 
tively. Integrating Eq. ([^ over all possible angles x ac¬ 
cepted by the detector, while averaging over the muon 
stopping volume results in a considerable reduction of 
the observed QI effect. Again, we create pseudo data for 


the real geometry, fit the resonances with a simple sum 
of Lorentzians, and determine the resulting shift 5qj of 
the line centers. We notice that taking into account the 
scattered photons at yf 90° and also the inhomogeneous 
muon stop probability, hardly affects the final results. As 
can be seen in Table [Hi these shifts are much lower than 
the experimental accuracy of a few % of the linewidth for 
all muonic atoms considered here [3I1IMI- 

IV. CONCLUSION 

We quantified the line shift caused by quantum inter¬ 
ference for /i^H, and /i^He’*' resonances, assuming 
first a pointlike detector. For /r^H, the resulting shifts 
are small. Hence, quantum interference cannot be the 
source of the proton radius puzzle, which requires a shift 
of the resonance in /i^H by four linewidths [SHE!]. On 
the other hand, the influence of quantum interference for 
some resonances of /r^H can be as large as 12 % of the 
linewidth for a pointlike detector. 

However, we verified that even for those large QI shifts, 
obtained assuming a pointlike detector, the angular aver¬ 
aging caused by the large acceptance angle of the photon 
detector and the size of the muon stop volume in the 
CREMA experiment significantly reduces this effect to 
negligible values at the present level of accuracy. 
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FIG. 4. (Color online) Simplified view of the CREMA target 
volume. Muons enter the gas target from the left (blue ar¬ 
row) . The gray cuboid indicates the muon stop volume which 
is illuminated by a laser pulse propagating inside a multipass 
cavity (red arrows). The incoming direction and polarization 
are indicated by fei and Si. The pointlike detector assump¬ 
tion treats the case where the muonic atom is at the center 
(red circle) and the scattered photon k-^ is emitted along 62 
and hits the top detector in the center (black circle). 
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Appendix A: DIPOLE MATRIX ELEMENTS 


The dipole matrix elements Dpl/Zji in Eq. (2) are eval¬ 
uated using standard angular reduction methods, which 
can start by simply expanding the product of the photon 
polarization and the position vector (e-y ■ r) in a spherical 


basis, i.e., 

= {/'F'rn'J'] e.y ■ r \/FmJ) = 

1 

^ {P'F'm'J'] rx \/FmJ) . 

A ^-1 

(Al) 

where /? contains all additional quantum numbers of the 
atomic state besides F, m and J. Following the geometry 
and nomenclatnre of Fig. 2, the spherical form of the 
(normalized) polarization vectors are given by 


(±1) _ -r ■ (cos Xi=l=^ sin xi) 


= T 


ATT 




(±1) _ (cos X 2 cos e±i sin X 2 ) -(0) _ 


= T 


"W 


, £2 = ~ *^08 X 2 sm I 


(A2) 


where x = Xi- 

The matrix elements of rx can be further simplified 
by making use of the Wigner-Eckart theorem [17] and 
considering the overall atomic state being the product 
coupling of the nucleus and electron angular momenta. 


i.e., 


\PFm) = ^ {JmjImi\Fm) |/m/) |/3Jmj) . (A3) 

rrijmj 

Here, the quantities {ji'mj-^j 2 'mj/\j^mj/) stand for the 
Clebsch-Gordan coefficients. After employing sum rules 
of Clebsch-Gordan coefficients we get 


{/'F'm'J' \rx\/FmJ) = {-if 1 ^ ) { } 

f Vtfi nrkfcif ir»n f t-. dV. lie onn q 1 i-rk rlr»oc nrk-l- rin flncj CTkiT 


where the notation [ji, j 2 , •■•] is equal to (2ji-|-l)(2j2 + l)--- • Since 
the reduced matrix element (J'||rA|| Ji) is given by 

(/3V'||r||/3J) = (-l)^'-i/2vTAJi(^^/2 J 

provided that L' + L + 1 is even, where L is the orbital angular momentum of the atomic state. Gombining Eqs. (A4) 

JP.JP p'. ^ ' 

and (A5) and rearranging the terms, the quantities Sfi,i and VL/j'' 


{/'J'\\r\m , (A4) 

rx does not act on the spin part of the wavefunction, 

'n'||r||n) , (A5) 


and 

= [J,^,F/^J[Ff,F„Jf,Ji] ^ Q - 1 / 2 ) ( 1/2 0 -1/2 ) { ji 

with 


P p ■ 

a) of Ecfs. (3) and (4) can be written 
Sf„i = (n/||r||ni.) {n/\r\\ni) = j r^RfR,,dr j r^R,,Ridr = ^ 


as 


(A6) 


/ 1 F,. 

I Jf 


F. 1 F, 

Ji 7 Jv 


qF ^ y- f Ff 1 F,. \ 1 F, \ 

^ ^ V —TTif X 2 rrii^ I \ —rriiy Ai rui I 

Ai,A2 Friiy \ J / \ / 


nF„ 

^FjF, ’ 

(A7) 

(AS) 


The functions R in Eq. (A6) stand for the radial non- relativistic wavefunctions, which for the case of 2s —>■ 
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TABLE III. Values of the coefficients ao, 02 and 62 corresponding to the parametrizations x) = 0,0 + a 2 P 2 (cos 7 ) and 

(^> X) = 62P2 (cos 7). 


I 

1/2 

F, 

0 

0 

1 

1 

1 

1 

F. 

1 

1 

1 

1 

1 

2 

J. 

1/2 

3/2 

1/2 

1/2 

3/2 

3/2 

F,. 

1 

1 

2 

2 

1 

3/2 

3/2 

3/2 

3/2 

3/2 

flo 

2/81 

4/81 

4/81 

4/81 

2/81 

10/81 

02 

0 

-2/81 

0 

0 

1/162 

-7/162 

&2 

-4/81 

-2/81 

-2/27 

-1/27 

1 

1/2 

1/2 

1/2 

2,12 

3/2 

4/729 

0 

-8/729 


1/2 

1/2 

3/2 

3/2 

1/2 

32/729 

0 

-32/729 


1/2 

3/2 

3/2 

- 

- 

40/729 

-4/729 

- 


3/2 

1/2 

1/2 

3/2 

3/2 

32/729 

0 

-32/3645 


3/2 

1/2 

1/2 

5/2 

3/2 

32/729 

0 

-32/405 


3/2 

1/2 

3/2 

3/2 

1/2 

4/729 

0 

-4/729 


3/2 

1/2 

3/2 

3/2 

3/2 

4/729 

0 

16/3645 


3/2 

1/2 

3/2 

5/2 

3/2 

4/729 

0 

-4^05 


3/2 

3/2 

1/2 

3/2 

3/2 

40/729 

0 

-128/3645 


3/2 

3/2 

1/2 

5/2 

3/2 

40/729 

0 

-28/405 


3/2 

3/2 

3/2 

5/2 

3/2 

32/729 

64/18225 

-112/2025 


3/2 

5/2 

3/2 

- 

- 

4/27 

-28/675 

- 


side of Eq. (|A6|) . The quantity is the ratio between 


2p —>■ Is gives the numerical result shown on the right previous angles by cosy = cosy sin 0 . Aj'j‘'( 0 ,x) and 

(^’ x) are parametrized in terms of this angle 7 
by oo + a 2 P 2 {cosj) and l> 2 -P 2 (cos 7 ), resp ecti vely. The 
respective coefficients calculated using Eq. (A7) are listed 
in Table Hill 


the muon-nucleus reduced mass and the electron mass. 

In case of incident linear polarized photons, the dipole 
radiation pattern of the scattered photon depends only 
on the angle 7 between the incident polarization and the 
direction of the scattered photon, which is related to the 
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